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Abstract: We use a diffusion galactic model to analyze the end of the Galactic cosmic ray spectrum 
and its mixing with the extragalactic cosmic ray flux. We analyze the transition between Galactic and 
extragalactic components using two different extragalactic models. We compare the sum of the diffusive 
galactic spectrum and extragalactic spectrum with the available experimental data. 
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Introduction 

The cosmic ray energy spectrum extends for many 
orders of magnitude with a power law index « 2.7. 
Along this range of energies, the three spectral fea- 
tures are known: the first knee at E w 3 PeV, 
the second knee at E « 0.5 EeV and the ankle, 
a dip extending from the second knee to beyond 
10 EeV. The nature of the second knee and of the 
ankle are still uncertain [1]; a possible interpreta- 
tion of the two features is the transition between 
the galactic and extragalactic components. At en- 
ergies between 10 17 — 10 18 eVthe galactic super- 
nova remnants (SNRs) are expected to become in- 
efficient as accelerators. This fact, combined with 
magnetic deconfinement should mark the end of 
the Galactic component of cosmic rays, although 
the picture could be confused by the existence 
of additional Galactic accelerators at higher en- 
ergy. On the other hand, at energies above the sec- 
ond knee, extragalactic particles are able to travel 
from the nearest extragalactic sources in less than 
a Hubble time. Consequently, the spectrum may 
present above 10 17,5 eV a growing extragalactic 
component that becomes dominant above 10 19 eV. 
The region between the second knee and the ankle 
could be the transition region between the galactic 
and extragalactic components. 
In this work we analyze the transition region com- 
paring the diffusive galactic spectrum from SNRs 
with two different models of extragalactic spec- 



trum, one in which only protons [2] are injected 
at the sources and another in which a mixed com- 
position containing heavy nuclei [3] is injected. 

Diffusion Galactic model 

We used the numerical diffusive propagation code 
GALPROP [4, 5] to reproduce the galactic spec- 
trum from SuperNova Remnants (SNRs). The dif- 
fusive model is axisymmetric. The propagation 
region is, in cylindrical coordinates, bounded by 
R = Bf, = 30 kpc and z = Zh = 4 kpc, be- 
yond which free escape is assumed. The propaga- 
tion equation is: 
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where ^(r, p, t) is the density per unit of total par- 
ticle momentum, q(r,p) is the source term, D xx 
is the spatial diffusion coefficient, p = dp /dt is 
the momentum loss rate and Tf and r r are the time 
scale of fragmentation and the time scale of ra- 
dioactive decay respectively. The diffusion coeffi- 
cient is taken as /3Dq(p/ pr>) 5 , where p is the par- 
ticle rigidity, Do is the diffusion coefficient at a ref- 
erence rigidity po and S = 0.6. The distribution of 
cosmic rays sources used is that of Galactic SNRs 
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Figure I: Diffusive Galactic spectrum, <I>g, 
with renormalization of the CNO group (Z:6-12), 
$cno and Q f me heayy component (Z : 19 - 26), 
<&q . Agreement with experimental data is obtained 
for k f land k ^0.8 -1.2. 



(in agreement with EGRET gamma-ray data) [4]. 
The injection spectrum is a a power law function in 
rigidity with a break at rigidity po, beyond which 
it falls exponentially with a rigidity scale p c : 
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where a = 2.05, p = 1-8 ^ and p c = 1.26 PF. 
Stable nuclei with Z < 26 are injected, with 
energy independent isotopic abundances derived 
from low energy CR measurements [5]. 

Detailed and realistic interstellar molecular (H 2 ), 
atomic (H) and ionized (HI) hydrogen distributions 
are used [6] . 

Diffusive Galactic spectrum 

The diffusive galactic spectrum has been normal- 
ized to match KASCADE data at - 3 x 10 6 
GeV. While with this renormalization our spec- 
trum agrees with JACEE and Sokol data at lower 
energies, beyond the knee the diffusive spectrum 
presents a strong deficit of flux. Since at E > 
10 7 GeV the composition is dominated by inter- 
mediate (Z : 6 - 12) and heavier (Z : 19 - 26) 
nuclei, we renormalize these components by a fac- 
tor of 2, which produces a good agreement with the 
experimental data (see Fig. 1). 

The renormalized diffusive galactic spectrum re- 
produces well the data up to E « 10 8 GeV, be- 



yond which the spectrum falls steeply because of 
end of the SNR acceleration efficiency. 

Extragalactic spectrum 

In order to study how the transition between the 
galactic and extragalactic components takes place, 
we compare the galactic spectrum originated in 
SNR with two different possible scenarios for the 
extragalactic component. 

In the first model [2], a pure proton extragalac- 
tic spectrum, accelerated by a homogeneous dis- 
tribution of cosmic sources, is considered. Local 
overdensities/deficits of UHECR sources affect the 
shape of the GZK modulation, but do not affect 
the low energy region where the matching with the 
Galactic spectrum occurs. Different cases of lo- 
cal overdensity/deficity of sources are considered 
within this model. 

In the second model [3], the extragalactic spec- 
trum is calculated for a mixed composition at injec- 
tion typical of low energy cosmic rays for different 
source evolution models in red shift. 

In both cases, the various parameters of the models 
are tuned to fit the available CR data at UHE and 
are, in that highest energy regime, experimentally 
indistinguishable at present. 

Combined spectrum: matching Galac- 
tic and extragalactic components 

In order to study how the transition between the 
Galactic and extragalactic components takes place, 
we subtract the combined theoretical (Galactic plus 
extragalactic) spectrum from the available data. 
Two different approaches are used. 

First, we try to match the experimental data by 
varying the normalization of the heavy Galactic 
component, while keeping constant the previous 
renormalization of the CNO group. The best re- 
produced spectrum for the two extragalactic mod- 
els are shown in Figs. 2, 3 and 4. In the case of 
the proton model, a discontinuity appears when the 
two spectra are added, regardless of the lower limit 
adopted for the extragalactic component: 10 8 GeV 
(Fig.2) or 5 x 10 7 GeV (Fig.3). The latter cor- 
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Figure 2: Galactive and Extragalactic spectrum 
matching for the proton models for a EG lower 
energy limit of 10 8 GeV. The sum of the renor- 
malized diffusive Galactic spectrum and different 
extragalactic spectrum models (&eg) is shown for 
different renormalizations of the heavy component 
The CNO group (^g^ ) of the diffusive 
Galactic spectrum ($g<) has been renormalized by 
a factor 2.2. Different cases of local overden- 
sity /deficit of sources are considered [2]: (1) uni- 
versal spectrum, (2) and (3) overdensity of sources, 
(4) and (5) deficit of sources . 

responds to cosmic accelerators operating for the 
entire Hubble time. 

For both, proton and mixed-composition models, 
there is a flux deficit above 10 8 GeV. The problem 
is much stronger for the mixed-composition model 
where, regardless of the luminosity evolution of 
EG CR sources, the total spectrum presents a large 
deficit of flux between 10 8 and « 3 x 10 9 GeV 
(Fig.4). 

In order to solve this flux deficit, the only way 
out seems to be the introduction of an additional 
Galactic component. We estimate this component 




Figure 3: Idem to Fig. 2 but for the proton model 
with a lower energy limit of 5 * 10 T GeV. 



Figure 4: Idem to Fig. 3 but for the Galactive and 
Extragalactic spectrum for the mixed-composition 
model. The CNO group (<$>% NO ) of the diffusive 
Galactic spectrum ($g) has been renormalized by 
a factor 2. Three source evolution models are con- 
sidered [3]: (a) strong, (b) SFR and (c) uniform. 

by subtracting the sum of the diffusive Galactic 
and extragalactic fluxes from a smooth fit to the 
world data. This method confirms us the need of 
the renormalization of the CNO group by a fac- 
tor « 2.2 and « 2 for the proton and mixed- 
composition models, respectively. In the case of 
the proton models (Figs. 5, 6), the observed deficit 
can be resolved with an additional heavy compo- 
nent, while in the mixed-composition models this 
is not enough and we need one more additional 
heavy component (Fig.7). The additional compo- 
nent common to both families of models is ob- 
tained with a shift in energy of a factor w 1.5 of 
the diffusive galactic heavy component, renormal- 
ized by a factor « 0.04 in the case of the mixed- 
composition models and < 0.03 for the proton 
models respectively. The second additional com- 
ponent is obtained in an analogous way but with an 
energy- shift factor of « 1.8 and a renormalization 
by a factor 5 x 10 -4 . The corresponding spectra 
are shown in Figs. 5, 6 and 7. 

Discussion 

We have analyzed the matching conditions of the 
Galactic and extragalactic components of cosmic 
rays along the second knee and the ankle. 

It seems clear that an acceptable matching of 
the Galactic and extragalactic fluxes can only be 
achieved if the Galaxy has additional accelerators, 
besides the regular SNR, operating in the interstel- 
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Figure 5: Galactive and Extragalactic spectrum 
for the proton model for an EG lower limit of 
10 8 GeV\ the additional heavy component (j>Q A 
and the total heavy component (<pQ + (f>Q A ) are 
also shown. 




Figure 7: Galactive and Extragalactic spectrum for 
the mixed-composition model: the two additional 
heavy components (/)q A , (^ga2 anc * me total heavy 
component (4>q + <j)Q A + 0^ 2 ) are snown - 

acceleration requirements imposed on the Galaxy 
and probably rapidly spinning inductors, like neu- 
tron stars, could be invoked to fill in the gap. If this 
were the case, it is very likely that photon emis- 
sion at TeV energies should uncover the sources. 
Changes in propagation regime inside the Galaxy 
at these high energies should manifest as a dipolar 
anisotropy if enough statistics were available. 



Figure 6: Idem to Fig. 5 but for an EG lower limit 
of 5 x 10 7 GeV: the additional heavy component 
<Pq A and the total heavy component (<pQ + (/)q A ) 
are also shown. 

lar medium. In the particular case of the proton 
model, only one additional component is required, 
and this could well represent the contribution from 
compact and highly magnetized SNR, like those 
occurring in the central, high density regions of 
the Galactic bulge or the dense cores of molecular 
clouds. It must be noted, however, that a perfectly 
smooth match seems unrealistic and that some dis- 
continuity, whose magnitude depends mainly on 
the time depth from which EG CR are able to ar- 
rive at our Galaxy, should be eventually observable 
in the combined spectrum. 

The matching of the mixed-composition model is 
more complicated. The Galactic spectrum has to 
be extended up to the middle of the dip and this re- 
quires, besides the previous additional component, 
another high energy Galactic component. The ori- 
gin of these cosmic rays pushes even further the 
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